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atlon  ( or  regression  ) from  tracking  occurs.  If, 
however,  the  input  is  not  changing  too  rapidly, 
the  hunan  will  reacquire  the  target,  and  continue 
tracking.  A study  of  this  type  of  adaptability 
is  considered  here  using  the  phase  plane  with  a 
statistical  analysis  performed  over  a family  of 
four  input  forcing  functions. 
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non-s tat ionary  effects  will  occur  in  the  closed 
loop  interaction  between  the  man  and  the  machine. 
Figure  (3)  illustrates  the  velocity  and  acceler- 
ation profiles  of  the  forcing  function  in  the 
azimuth  axis  for  the  first  input  considered  here. 
It  is  observed  that  large  variations  of  the  velo- 
city and  acceleration  profiles  may  give  rise  to 
adaptive  behavior  of  the  man  in  the  loop  as  he  is 
required  to  perform  this  type  of  tracking  task. 
Figure  (4)  Illustrates  a segmentation  of  the 
tracking  task  which  has  been  observed  to  occur 
with  humans  as  they  are  forced  to  track  this  type 
of  input.  From  the  time  period  tQ  to  t^ , the  vel- 
ocity and  acceleration  of  the  input  are  sufficien- 
tly small  in  magnitude  such  that  the  tracker  can 
maintain  the  plant's  output  to  follow  the  input 
forcing  function.  From  the  time  period  t^  to  t2> 
the  input  is  too  difficult  to  follow  and  the  hu- 
man regresses  and  no  longer  attempts  to  follow  the 
input.  During  this  period  the  error  signal  in 
the  closed  loop  system  is  dominated  by  the  signal 
f(t)  and  the  system  is  essentially  open  loop.  At 
the  time  t2  the  velocity  and  acceleration  of  the 
input  are  reduced  in  magnitude  such  that  the  human 
can  then  reacquire  the  target. 

The  experimental  data  used  here  came  from  a 
simulation  study  of  male  and  female  trackers 
which  consisted  of  5 teams  of  2 trackers  each  who 
were  trained  for  3 weeks  prior  to  the  data  presen- 
ted here.  The  training  period  consisted  of  24  runs 
a day  on  4 different  simulated  flight  trajectories 
(input  forcing  functions).  The  deterministic  tra- 
jectories were  presented  to  the  subjects  on  a 
random  basis.  During  the  experimental  period  the 
subjects  tracked  12  runs  per  day.  The  four  forcing 
functions  used  in  this  study  resembled  the  shape 
of  figure  (2)  with  differences  in  the  shape  and 
magnitudes  of  the  velocity  and  acceleration  pro- 
files. These  forcing  functions  were  chosen  based 
on  known  aircraft  maneuvers  which  were  of  interest 
The  plant  dynamics  which  appear  in  figure  (1)  can 
be  represented  by  the  following  lumped  transfer 
function: 
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1.  Introduction 


The  study  of  a human  Involved  in  a closed 
loop  tracking  situation  in  which  the  task  is  to 
control  the  dynamical  response  of  some  mechanical 
device  has  been  a problem  of  interest  for  many 
years.  The  use  of  phase  plane  analysis  to  study 
man-machine  interaction  is  not  nev  A has  been 
considered  previously,  for  example  oy  Phatak  and 
BekeyDJ  . The  approach  considered  here  differs 
from  previous  approaches  because  here  the  problem 
centers  on  the  identification  of  the  human  limita- 
tions and  also  on  the  performance  index  eesodated 
with  the  man-machine  interaction.  Figure  (1)  ill- 
ustrates the  typical  man  in  the  loop  problem  con- 
sidered in  this  paper  for  the  two  operators  in- 
volved in  e two  dimensional  tracking  task.  Figure 
(2)  illustrates  one  of  the  four  input  forcing 
functions  in  the  horizontal  and  vertical  axis  for 
this  two  dimensional  task.  Each  operator  was  re- 
quired to  track  the  forcing  function  in  hie  axis 
(azimuth  or  elevation  which  correspond  to  horizon- 
tal and  vertical,  respectively).  This  type  of  sim- 
ulation is  typical  of  AAA  (anti-aircraft-artil- 
lery) simulations  as  considered  in  [2 ,3]  . As  a 

result  of  this  type  offorclng  function  many 
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In  order  to  study  these  tracking  problems  with 
this  non-s tationary  behavior,  it  is  of  interes 
to  observe  these  effects  on  behavior  in  the 
phase  plane. 


2.  A Phase  Plana  Analysis  of  The  Tracking  Task 


Figure  (5)  Illustrates  phase  plane  trajectories 
of  the  closed  loop  error  signal  during  typical  seg- 
ments of  tracking  (from  tQ  to  tj) , during  the  re- 
gression period  (from  t^  to  t2),  and  during  the 
reacquisition  of  the  target  and  final  tracking  per- 
iod (from  t2  to  tf) . During  the  time  period  the 
human  is  seriously  tracking  the  error  signal,  the 
phase  plane  is  elliptical  in  shape.  During  the  re- 
gression and  reacquisltion  period  the  trajectories 
spiral  out  from  the  origin  (and  return  to  the  ori- 
gin) as  if  the  system  were  first  unstable  and  then 
stable.  Tracking  behavior  in  the  final  segment  is 
elliptical  in  shape  similar  to  the  behavior  in  the 
initial  tracking  segment.  The  elliptical  shape  in 
the  phase  plane  has  a physical  Interpretation  which 
implies  that  greater  importance  is  placed  on  re- 
ducing the  position  value  of  the  error  signal  with 
less  importance  placed  on  reducing  the  derivative 
of  the  error  signal.  One  may  interpret  figure  (5) 
in  a performance  index  sense.  Consider  the  follow- 
ing performance  index: 

J - /T  [e2 ( t)  % + 62(t)  R^dt  (2) 

with  R]>0,  R2>0,  and  T is  the  period  of  time  it 
takes  an  ellipse  to  encircle  the  origin.  It  is 
noted  that  the  scalar  variables  e(t)  and  e(t)  are 
weighted  but  no  penality  weighting  is  assigned  to 
the  control  vector  u(t)  generated  by  the  man.  The 
performance  index  (2)  implicitly  weights  the  con- 
trol vector  u(t)  due  to  the  closed  loop  dependence 
between  e(t)  and  u(t).  The  elliptical  shapes  which 
appear  in  figure  (5)  imply  that  R,>  Rj  or  that  more 
penality  is  associated  with  minimizing  the  position 
state  of  the  error  signal  and  less  penality  is 
associated  with  minimizing  the  rate  of  the  error 
signal.  An  approach  will  now  be  introduced  to  in- 
vestigate the  weightings  Rj  and  R2  in  equation  (2) 
and  to  test  statistical  hypothesis  concerning  them. 

3.  An  Approach  To  Identify  The  Performance  Index 


In  order  to  identify  a performance  index  of  the 
form  (2) , consider  the  equations  of  an  ellipse  with 
respect  to  the  variables  e(t)  and  e(t): 


e2(t)  e2(t) 

R2  *1 


1 (3) 


or  e2(t)  Ri  f e2(t)  R2  “ R^  R2  (4) 


But: 


Constant 


R2 


dt 


[e2(t)  R-l  + e2(t)  R2] 


dt 


(5) 

J (6) 


Therefore  J can  be  determined  knowing  only  R^,  R2, 
and  T.  Since  J is  explicitly  expressed  by  equation 
(6)  over  the  time  period  [O.Tj,  then  it  is  required 
to  identify  T as  the  time  necessary  for  the  ellipse 
to  encircle  the  origin.  Therefore  it  is  only  nec- 
essary to  identify  the  variables  R^,  R2t  and  T dur- 
ing the  various  stages  of  tracking  and  also  to  study 
how  these  variables  change  as  a function  of  the 


velocity  and  acceleration  characteristics  of  the 
input  forcing  functions.  The  consistency  (measured 
by  means  and  variances)  of  the  performance  index 
parameters  R2,  R2,  and  T can  be  investigated  as  a 
function  of  different  velocities  and  accelerations 
of  the  4 different  input  forcing  functions.  From 
equation  (3)  , R^  and  R2  become  the  axis  of  the 
ellipse  and  the  time  period  T can  be  read  off  para- 
metrically in  the  phase  plane.  Regression  was  de- 
fined as  a trajectory  whose  distance  from  the 
origin  is  in  excess  of  twice  the  largest  axis  of 
ellipse. 

It  is  interesting  to  investigate  if  the  regress- 
ion from  tracking  is  due  to  the  velocity  or  accel- 
eration of  the  input  or  possibly  to  both  effects 
taken  together.  The  following  statistical  section 
will  describe  the  manner  in  which  these  type  of 
hypothesis  can  be  investigated. 

4.  A Statistical  Analysis  of  The  Adaptation  Effects 

In  an  effort  to  identify  the  performance  index 
of  the  form  (2)  from  the  data,  it  is  necessary  from 
equation  (4)  to  identify  Rj^,  R2,  and  T from  the 
elliptical  phase  plane  plots.  The  variables  R^,  R2, 
and  T will  vary  over  the  initial  and  final  stages 
of  tracking  and  also  over  all  four  forcing  functions. 
The  following  table  is  constructed  based  on  data 
from  as  many  as  18  replications  of  the  four  forcing 
functions  for  the  best  team  in  this  study.  The 
criteria  for  the  selection  of  the  best  team  from 
the  five  possible  candidate  teams  was  based  on  the 
maximum  time  during  the  45  second  run  that  the 
error  signal  was  within  a specified  error  window 
size.  Table  (1)  illustrates  the  results  obtained 
for  this  experiment  during  the  time  perlod£t0  , tjJ . 


Table  (1) -Performance  Index  Coefficients  £t0  , t J 


FF#1 

| FF#2 

' FF#3 

| FF#4 

ca 

MU 

o 

wum 

d 

IMiM 

C9 

IM21 

R^mean) 

sta 

gjj 

U 

U 

111 

m 

HI 

Rl(e.d.) 

.02 

.01 

.07 

.02 

.04 

.03 

.04 

.02 

R2(mean) 

.06 

.03 

.12 

.04 

.06 

.03 

.07 

.03 

R2(s.d.) 

.01 

.01 

.03 

.01 

.03 

.01 

.02 

.01 

T(mean) 

m 

2.5 

3.3 

2.8 

3.2 

2.9 

3.2 

2.8 

T(s.d.) 

.93  | 

.48 

.88 

.50 

IQ 

.70 

El 

.75 

where  Az  and  El  stand  for  Azimuth  and  Elevation  axis, 
respectively .The  same  results  are  determined  for 
the  time  period  £t2  > t^J  : 
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Table  (2) -Performance  Index  Coefficients  £t2  , tj] 


R^(mean) 


^(s.d.) 


R2 (mean) 


1^(8. d.) 


T(mean) 


T(a.d.) 


FF#1  FF#2  I FF#3 


.29  1 .11  1 .29  .10  .23  .14  .30  I .11 


.03  .03  .11  .02  .04  .03  .04  .02 


.07  .03  .10  .03  .07  .03  .07  .03 


.02  .01  .04  .01  .04  .01  .02  .01 


3.2  2.6  3.3  2.8  3.2 


.64  I .68  I .96  .50  1.1  .57  .69  I .72 


Table  (5)  was  computed  from  the  data  in  order 
to  study  velocity  and  acceleration  effects  in  the 
reacquisition  of  the  target: 

Table  (5)  - An  Investigation  of  Reacquisition 
Characteristics  Over  The  4 Inputs: 


FF#1 

FF#2 

FF#3 

FF#4 

t2 (mean)  None  22.0  30.1  28.1  33.2  30.0  30.0  27.5 


t„(s.d.)|  None  I 1.5  I 2.2 


vel.  at 
t2(mean) 


acc.  at 
t2(mean) 


-.60  -3. 


-3.0  -3.5 


-.23  1.16  -.13  .2  | -.20  0.7  .18 


It  is  necessary  to  now  teat  whether  the  track- 
ing performance  index  that  was  used  in  the  closed 
loop  system  changed  significantly  from  the  time 
period[t0  , tH  to  [t2  , tjj  . Table  (3)  represents 
this  aspect  of  changes  in  tracking  strategy  of  the 
human: 


Table  (3)  - An  Investigation  of  Changes  in  The 
Performance  Index  prior  to  and  after  Regression: 
( Fly-By  # 3) 


Mean  S.  0.  Mean 
Prior  Prior  After 
to  ti  to  ti  t? 


.1318  .0275  .1355  .0266 


.0356 


stat.  I level 


Finally  to  study  the  variability  of  velocity  and/or 
acceleration  as  the  primary  factor  in  causing 
regression,  t-tests  were  made  against  zero  for  the 
mean  (over  all  four  forcing  functions)  of  the  velo- 
cities versus  zero  and  the  accelerations  versus 
zero.  These  t-tests  were  also  conducted  for  the 
reacquisition  task  and  the  results  are  displayed 
in  table  (6).  A higher k level  indicates  that  this 
variable  has  greater  variance  if  the  human  it 
using  this  variable  in  his  decision  making  process. 

In  other  words,  if  the  human  is  a minimum  variance 
estimator,  he  will  use  the  variable  which  gives 
rise  to  the  smallest  «4  level. 

Table  (6)  - An  Investigation  of  Velocity  and 
Acceleration  Effects  At  Regression  and  Reacuulsltion: 


t <4 

stat.  level 


2.02  .114 


.0277  .0085  .0309  .0089 


-.966  1.343 


El  2.88  .70  2.78  0.57  14  0.4145  .682 


In  order  to  determine  the  source  of  regression  (vel 
ocity  and/or  acceleration),  table  (4)  was  computed 
from  the  data  to  study  this  effect: 

Table  (4)-  An  Investigation  of  Regression  Character 
is tics  Over  The  4 Inputs: 


1.83  .116 


2.48  .068 


8.80  . 0001 


t2(mean) 

Az 

None 

El 

17.5 

Az 

11.8 

El 

14.1 

Az 

14.5 

El 

12,5 

Az 

14.0 

El 

loTT 

tj.(s.d.) 

None 

0.5 

2.47 

1.1 

2.40 

1.47 

1.2 

2.54 

vel.  at 
tx(mean) 

1 - 

.68 

' 

.40 

.50 

-1.3 

.25 

-3.1 

2.75 

see.  at 
tx(mean) 

- 

-.23 

-.70 

-.10 

-.90 

i 

-.60 

-.60 

-0.8 

Summary  and  Conclusions: 


The  results  of  tables  (3),  (5),  and  (6)  indicate 
the  following  (based  on  the  «<  levels  and  the  data 
base  used  here) : 

(1)  The  performance  index  which  describes  tracking 
prior  to  and  after  regression  does  not  change 
statistically  in  its  weighting  coefficients  or 
the  time  necessary  for  the  ellipse  to  encircle 
the  origin.  In  other  words,  the  tracking 
prior  to  regression  and  after  reacquisition 
remains  the  same  in  a statistical  sense. 
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(2)  There  Is  a significant  less  variability  if  the 
human  were  to  use  acceleration  of  the  target  as 
his  decision  mechanism  in  the  regression  from 
the  task  and  reacquisition  of  the  task  as  com- 
pared to  the  velocity  of  the  target.  In  other 
words, the  factors  which  determine  regression 
and  reacquisition  seem  to  be  more  closely 
tied  to  acceleration  (with  less  variance)  as 
compared  to  the  velocity  of  the  target. 
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